Observation of an exceptional point in a chaotic optical microcavity 
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We present spectroscopic observation of an exceptional point or the transition point between dia- 
batic crossing and avoided crossing of neighboring quasi-eigenmodes in a chaotic optical microcavity 
with a large size parameter. The transition to the avoided crossing was impeded until the degree 
of deformation exceeded a threshold deformation owing to the system's openness also enhanced by 
the shape deformation. As a result, a singular topology was observed around the exceptional point 
on the eigenfrequency surfaces, resulting in fundamental inconsistency in mode labeling. 

PACS numbers: 05.45.Mt,42.55.Sa,42.65.Sf 



Dielectric optical microcavities are widely used in var- 
ious optoelectronics applications such as add/drop fil- 
ters, low-threshold microlasers [I], single- molecule sen- 
sors tunable optical frequency comb 0] and optome- 
chanical oscillators owing to the high quality factor Q 
of supported modes. In addition, asymmetric optical mi- 
crocavities have drawn much attention because they ex- 
hibit directional output emission of high Q modes. They 
can also serve as a useful platform for investigating the 
correspondence between quasi-eigenstates and associated 
chaotic classical dynamics in mesoscopic systems due to 
the well-known one-to-one correspondence between the 
Schrodinger equation and the Maxwell equation in bil- 
liard problems @, 0, S, 0, Efll, EH, El, El ■ Level dy- 
namics of interacting modes have also been studied for 
tailoring output directionality and quality factors of as- 
sociated modes [LJ, [l5( in microcavities. 

One of the key issues related to the level dynamics 
is an exceptional point (EP), a singular point generally 
existing in a parameter space of a non-Hermitian sys- 



tem [lfl [17J, |18j. EPs have been observed in various 
systems such as acoustic systems atoms in optical 
lattices [2(J, and complex atoms in laser fields 21 1. In 



coupled microwave cavities the EP was studied in terms 
of resonance mode distributions H, 23| . Existence of an 
EP has been predicted for mesoscopic systems such as 
a Rydberg atom in a strong magnetic field 24[ and in 
a stadium-shape microcavity [25| . where classical chaos 
is also important. Yet their experimental verifications in 
these classically chaotic systems have not been reported. 

In a recent work of Ref. [26[, we observed both di- 
abatic crossing (DC) and avoided crossing (AC) in a 
chaotic optical microcavity (COM), depending on the 
degree of cavity-shape asymmetry. However, it was not 
possible to observe an EP, where DC and AC coalesce in 
a two-dimensional parameter space. The problem was a 
finite spectral resolution, which prevented us from telling 
whether two modes were crossing each other or they were 



avoiding each other with a gap smaller than the spectral 
resolution. Similar problems were also found with in- 
distinguishable static envelops of diabatic and avoided 
crossing modes in a high-Q toroidal microcavity 27|. 

In this Letter we report the first observation of an EP 
in a high-Q asymmetric microcavity or COM with its 
boundary shape continuously variable. Our experiment 
done in a single COM, not in coupled cavities as before 
[22L [H, was made possible by introducing an internal 
parameter, a quasi-continuous variable in a semiclassi- 
cal regime of large size parameter, and is based on the 
fact that two modes undergoing an AC exhibit funda- 
mentally different output coupling signatures from those 
of DC modes. Moreover, we elucidate the resulting sin- 
gular topology in the eigenenergy surfaces around the EP 
as an inherent source of impossibility of consistent mode 
labeling in open mesoscopic systems. 

Our experiment was performed with a COM or a 
two-dimensional microcavity formed by a liquid jet of 
ethanol (refractive index m=1.361 at 610 nm) dilutely 
doped with Rhodaminc dye molecules. The boundary 
profile of the microcavity is approximated by r(4>) ~ 
a(l + rj cos 2(f) + ei] 2 cos 4<j>) in the polar coordinates with 
a ~14.9±0.1/mi and e = 0.42 ± 0.05 [HIU. The defor- 
mation parameter 77 can be continuously varied from 0% 
to 26% at will. The size parameter, defined as 27rma/A 
with A the wavelength, is about 190 for A=660 nm. 

The COM was pumped by a cw argon-ion laser and 
the cavity-modified fluorescence (CMF) or lasing light 
from the COM was measured with a spectrometer. The 
polarization direction for both pumping and detection 
was parallel to the COM column 0, EH, 0]- There were 
observed five different mode sequences recurring with an 
interval or a free spectral range (FSR) of about 2.3 THz. 
Mode label I is assigned to each mode sequence by the 
scheme explained in Ref. [26| and the recurring modes in 
a mode sequence are indexed by mode number n. 

In Fig. m we observe 1=2 and 1=4 modes undergoing 



2 



Internal parameter n 

181 




465 468 

Frequency (THz) 



1 1 1 1— 

» 


— I 1 1 1 ' 

+ ■ 




y - 


. ' 1 1 





16 18 20 22 

Deformation parameter, ti (%) 




10 12 14 16 18 20 22 

Deformation parameter, r) (%) 



FIG. 2: (a) Magnitude Aza4 of AC between two uncoupled 
states, l§ = \ and £o=4 modes, as a function of r\. Solid line is 
a linear fit. (b) Decay rates 71 and 74 of Zo=l,4 uncoupled 
states in the spectral region where AC or DC occurs. Dotted 
line for 71 is an extrapolation beyond spectral visibility. 



FIG. 1: Eigenfrequencies of 1=2 (in orange) and £=4 (in blue) 
modes as a function of the deformation rj and the internal 
parameter n. The magnitude of AC increases as r\ increases. 
The magnified view shows the high resolution spectra taken 
with a spectrometer of 0.05 nm resolution. Green curves are 
Lorentzian fits. 

ACs as a function of 77 at ra=179, 180, 181 and 182, re- 
spectively, indicating that the separation of two modes 
can be adjusted continuously by changing 77 for a fixed n. 
We also observe ACs as we move horizontally in Fig.Q]for 
a fixed 77, for example 77=0.176, from one FSR to another. 
This is equivalent to scanning the mode number n, a dis- 
crete internal parameter. The relative frequency between 
1=2 and 4 modes is changed by 0.083 THz or 3.6% (when 
they are well apart) of their FSRs when n is shifted by 
one. From the locations of ACs, (n, ry)=(179, 0.162), 
(180, 0.167), (181, 0.172), (182, 0.176), we also find that 
the same 0.083 THz induced by n shift can be continu- 
ously scanned by changing 77 by about 0.005. Therefore, 
by coarse scanning of n shift for a fixed 77, followed by a 
fine adjustment of 77 within 0.005, we can bring any two 
modes together, far beyond our spectrometer resolution 
(^50 GHz). Exact definition of the internal parameter n 
is given below when we discuss Fig. |H 

Before we examine the interaction between modes 1=2 
and 4 closely, let us recapitulate a theoretical background 
on mode-mode interaction. It was shown in Ref. [26j that 
the mode-mode interaction in a COM can be described 
by a 2-by-2 non-Hcrmitian symmetric Hamiltonian. Its 
diagonal elements, given by Ej(n,rf) = i/j(n,rj) — 77^(77) 
(j=a, b) (with the Planck constant h = 1), are the eigen- 
values of quasi-eigenmodes when their eigenvalues are 
well separated, i.e., effectively uncoupled. The imagi- 
nary part corresponds to the decay rate of the mode. 
The symmetric off-diagonal element, denoted by C(rj), is 
the mode-mode coupling constant induced by the nonin- 
tegrability or the cavity shape asymmetry while its de- 
pendence on the internal parameter n is negligible in the 
small spectral range of interest. 

The uncoupled quasi-eigenmodes of different radial 
mode order have different FSRs. Thus, by shifting the 



internal parameter n for a fixed 77 followed by a fine 
adjustment of 77 as explained above, we can bring any 
two uncoupled modes together and make the mode-mode 
coupling come into play. In this case the system has 
new eigenvalues E± given by E± = (E a + E b )/2 ± 
^(E a - E b ) 2 /4 + C 2 . When v a = v b {= u ), the new 
eigenvalues are given by E± = v§ ± (C 2 — 7 2 .) 1 / 2 — 77+ 
with 7± = |7 a ± 7&I/2. When C < 7_, a splitting occurs 
in the imaginary part of energy, corresponding to a cross- 
ing in the real part as we vary n. When C > 7_, an AC, 
a splitting in the real part, would occur as n is varied. 
This splitting is analogous to the normal mode splitting 
in coupled cavities [22| or in atom-cavity systems [30j ]. 
However, in our COM the internal coupling C is induced 
by the cavity shape asymmetry, related to chaotic ray 
transport between phase-space regions associated with 
the involved quasi-eigenmodes [iij ]. 

It should be noted that 1=2 and 1=4 modes in Fig. [1] 
correspond to the new quasi-eigenvalues E±, the result 
of AC between uncoupled j=a, b quasi-eigenmodes. We 
have systematically measured the sizes of AC between 
1=2 and l=A modes for various cavity deformation rang- 
ing from 0.14 to 0.23. We can also obtain the decay 
rates or the half linewidths of the uncoupled states in 
the spectral region where AC or DC occurs. In fact, we 
have identified from the spectrum evolution that these 
uncoupled states have evolved from whispering gallery 
modes (WGMs) with radial mode order Z =l and 4 in 
a circular cavity as 77 is gradually increased (26[. Thus 
we can label these uncoupled states j=a,b above by the 
same Zo=l and 4 as those of the original WGMs. The 
measured magnitude Az/14 of AC between lo=l and 4 
states and their decay rates 71 and 74 are summarized in 
Figs. [5]Ja) and[^b), respectively. 

It is interesting to note that the observed At^'s ap- 
pear to be well fit by a straight line with an 77-axis off- 
set, indicating there exist a threshold deformation for the 
splitting. The transition from DC to AC is suppressed 
up to a threshold deformation due to the openness which 
is also enhanced by system's nonintegrability. This fea- 
ture is quite interesting since the nonintegrability in an 
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open chaotic billiard induces both coupling and open- 
ness (the latter summed up by decay rates) to grow and 
it is not obvious how rapidly the coupling and the decay 
rates would increase as the degree of nonintegrability r\ 
rises. In our example of Zo— 1 and 4, the coupling (ini- 
tially zero) catches up the decay rates (initially finite) 
as X] grows. Contrarily, in a closed system an AC would 
occur as long as 77>0 however small it is since there is no 
decay by definition. 

From the 77-axis offset of a linear fit in Fig. [U(a) , one 
may expect that the two modes Iq— 1 and 4 would un- 
dergo a transition from AC to DC, or vice versa, at 77 ~ 
0.15. However, the linewidth (<3 GHz) of these modes 
become much narrower than our spectral resolution once 
r] is decreased below 0.16, and thus we cannot determine 
the actual value of threshold deformation by simply not- 
ing the disappearance of the splitting between the two 
modes. In fact, this problem would persist no matter 
how high the spectral resolution is unless it is infinite. 

One can solve this rather inherent problem by realiz- 
ing the fact that two modes undergoing a DC maintain 
their original linewidths whereas two modes undergoing 
an AC share a common linewidth given by a geometric 
mean of the original linewidths. This fundamental differ- 
ence is then translated to substantially different output 
coupling efficiencies for those two cases as to be seen be- 
low. Here the output coupling efficiency is a measure of 
the mode strength seen in the spectrum, given by the 
ratio £j=lj/{lj + 7abs) for the jth mode with 7 a b s the 
absorption rate of the cavity medium [3l[ . 

For Zo=l and 4 modes, we have (7i,74,7abs) ~ 
(2xl0- 6 ,4xl0- 3 ,2xl0- 5 ) THz and thus 71 < 7 abs < 
74 in the spectral region where the transition is expected, 
and therefore e\ ~ 11/ labs *Cl and €4 ~1. This indicates 
that Zo=l mode is hardly visible while /rj=4 mode has 
almost maximum visibility in the spectrum when these 
two modes are well separated. This behavior can be 
seen in Figs. EIa)-D3c), where Iq=1 modes are not visible 
in the fluorescence spectra whereas they become promi- 
nent when the pump power is increased beyond its lasing 
threshold as seen in Fig. [3fd). 

When an AC occurs (C > 7_ ~ 74/2), the two 
modes form a doublet with each having a decay rate of 
7+ ~ 74/2. The output coupling efficiency of each is then 
given by e + = 7+/ (7+ + labs) ~ 1- When the two modes 
barely undergo an AC, they appear to be overlapped into 
a single peak in the spectrum with its output coupling 
efficiency given by £ac = 2e+. This effect is clearly ob- 
served in Fig.[3ja), where the peak at 664.5 nm is twice 
larger than the other peaks (Zrj=4 modes) marked by blue 
arrows, indicating an AC takes place there between lo=l 
and 4 modes. 

When the two modes undergo a crossing (C < 7_), 
on the other hand, decay rate of each mode is modified 
from its original value as 7/4 = 7+ — (7^ — C 2 ) 1 ' 2 and 
7/4 = 7+ + (7 2 — C 2 ) 1 / 2 . The resulting output coupling 
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FIG. 3: CMF spectra for the case of (a) 77=0.125, (b) 77=0.102 
and (c) 77=0.078. Orange and blue arrows indicate Zo=l, 4 
modes, respectively. The output coupling efficiency of £rj=l 
modes is too small for these modes to be observed in (a)- 
(c). Red arrows indicate where /,o = l and 4 modes appear to 
be overlapped. Deformation r\ was fine tuned within 0.005 
to maximize the peak, (d) Pump power is increased above 
a lasing threshold for io=l mode at 661 nm with 77=0.120. 
2rj=l modes are clearly seen, (e) Observed output coupling 
efficiency (black squares) and the expected output coupling 
efficiency eac (red squares). Slight reduction of cac itself 
is mostly due to the decrease of decay rates of the involved 
modes. Black dotted line is a spline fit for visual guidance. 



efficiency when the two modes are overlapped is then 
given by e DC = 7 / i/(7 / i + labs) + lU/ilU + labs) < 
2e + , where the equality holds when C = 7_ or at the 
transition point. Particularly, if C <C 7- , 7'i ~ 71 
and 7/4 ~ 74, and the corresponding output coupling 
efficiency is approximately given by cdc Therefore, 
the transition from AC to DC is signaled by a substantial 
reduction of the output coupling efficiency from 2e + , as 
seen in Figs. G2h) and[3](c), when two modes appear to 
be overlapped. 

We examined the output coupling efficiency (black 
squares) for Zo=l and 4 modes when they appear to be 
overlapped. We observed that this coupling efficiency 
starts to deviate noticeably from the expected output 
coupling efficiency eac (red squares) as the deformation 
parameter is reduced below 0.125 as shown in Fig. |3je) . 
From this observation we conclude that the two modes 
undergo a transition from AC to DC, or vice versa, at 
77o=0.125±0.005. 

A parameter-space point at which the transition from 
DC to AC takes place or E± collapse to one is an EP, 
a topological singular point. The singular nature of the 
EP is revealed when we examine the eigenfrequency sur- 
faces of 1=2 and 4 modes in a 71-77 parameter space. The 
eigenfrequency surfaces E± (n, 77) are constructed in the 
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FIG. 4: (a) The observed spectrum for 77=0.187. (b) The rel- 
ative frequencies of observed quasi-eigenmodes with respect 
to the reference frequency of the same n as a function of n. 
(c) Eigenfrequency surfaces of two coupled states 1=2 (orange 
dots) and 1=4 modes (blue dots) in n-rj parameter space show 
a complex-square-root-like topology with a branch-point sin- 
gularity at (no, 7)o)=(175, 0.125). 



following way. We first define reference frequencies as the 
resonance frequencies of /o=3 WGM in a circular cavity 
whose round trip length is the same as that of the COM 
under investigation, as shown in Fig.|3Ja). The spectrum 
is then evenly divided into segments with each being as 
wide as one free spectral range of the Z =3 WGM so that 
each segment is indexed as n of that WGM. This n is pre- 
cisely the internal parameter used throughout this work. 
Its value, representing the number of wavelengths fitting 
the orbit associated with the mode, is as large as the size 
parameter. The relative frequencies of observed quasi- 
eigenmodes assigned with n, with respect to the reference 
frequency of the same n, are then plotted as a function 
of n as shown in Fig. H^b). By repeating this procedure 
for other rj and then combining the results as a function 
of n and r/, we finally obtain the eigenfrequency surfaces 
E±(n, rj) as shown in Fig.lUc). The resulting surfaces ex- 
hibit a complex-square-root-function-like topology with a 
branch-point singularity at the EP located approximately 
at (no, 7/- )=(175, 0.125), where both an AC and a DC be- 
tween Zo=l and 4 modes coalesce. 

The singular topology of the eigenenergy surfaces 
around the EP result in a fundamental inconsistency in 
assigning mode labels to quasi-eigenmodes. In order to il- 



lustrate this point, let us consider a cyclic variation of (n, 
rf) as shown in Fig. Egc): (170, 0.18)^(190, 0.18)^(190, 
0.10)^(170, 0.10)^(170, 0.18) enclosing the EP. If we 
choose 1=2 mode at (170, 0.18) and follow the mode un- 
der this cyclic variation, we end up with a different 1=4 
mode in the end after traversing A-^B^C^D^E on 
the energy surface. What happens is that the mode la- 
bel abruptly changes from 1=2 to 1=4 when we pass by 
the EP during adiabatic process of increasing -q (D^E), 
and moreover, there is no way to avoid this inconsis- 
tency however differently we label the modes. We have 
to perform the cyclic variation once more, traversing 
E^F^G— >H— >A around the EP, in order to come back 
to the same starting mode. This consideration reveals 
a fundamental inconsistency in assigning mode labels to 
quasi-eigenmodes in nonintcgrable open systems. This 
ambiguity is a direct consequence of the singular topol- 
ogy around the EP. 

In conclusion, we have observed an EP or the transi- 
tion point between DC and AC in a COM by utilizing 
different output coupling efficiencies in those two cases. 
The observed quasi-eigenfrequencies of interacting modes 
exhibit a branch-point topology, which is the very origin 
of impossibility of consistent mode labeling in these open 
chaotic systems. Our spectroscopic method using output 
coupling efficiency can be applied to other optical and 
optomechanical systems to resolve DC/ AC static mode 
envelops [27j ]. 
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